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■ ABSTRACT 

Of 

Q . We present the results of the Suzaku observation of the Seyfert 2 galaxy NGC 4507. 

■ This source is one of the X-ray brightest Compton-thin Seyfert 2s and a candidate for a 
I variable absorber Suzaku caught NGC 4507 in a highly absorbed state characterised by 

■ a high column density (iVn 8 x IQ^""* cm~^), a strong reflected component {R ^ 1.9) 

and a high equivalent width Fe Ka emission line {EW ~ 500 eV). The Fe Ka emission 
Une is unresolved at the resolution of the Suzaku CCDs (cr < 30 eV or FWHM < 3000 

^ ■ km s~^) and most likely originates in a distant absorber The Fe K/3 emission line is also 

Q>^ I clearly detected and its intensity is marginally higher than the theoretical value for low 

Tij- . ionisation Fe. A comparison with previous observations performed with XMM-Newton and 

CTS ■ BeppoSAX reveals that the X-ray spectral curvature changes on a timescale of a few months. 

\ We analysed all these historical observations, with standard models as well as with a most 

. recent model for a toroidal reprocessor and found that the main driver of the observed 2-10 

keV spectral variability is a change of the line-of-sight obscuration, varying from ~ 4 x 10^"^ 
psj , cm^^ to ~ 9 X 10^^ cm^^. The primary continuum is also variable, although its photon 

index does not appear to vary, while the Fe Ka line and reflection component are consistent 
with being constant across the observations. This suggests the presence of a rather constant 
reprocessor and that the observed line of sight A^h variability is either due to a certain degree 
of dumpiness of the putative torus or due to the presence of a second clumpy absorber. 
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1 INTRODUCTION 

The main ingred ient of the widely accepted Unified Model 
jAntonucci 1 1 19931) of Active Galactic Nuclei (AGN) is the pres- 
ence along the line of sight towards type 2 (or obscured) AGN of 
optically-thick material covering a wide solid angle. This absorber 
was thought to be rather uniform and distributed in a toroidal 
geomet ry, whic h is locate d at a pc scale distance from the nuclear 
region Jurrv & Pad ovanii 1 19951) or in the form of a bi-conical 
outflow ( lElvisll2000l) . Although this model has been to first order 
confirmed, it is now clear that this is only a simple scenario, 
which does not hold for all AG N jBianchi et al.ll201^ : [Turner et al.l 
I2OI2I : [Turner & MilleJ |2009|) . In particular X-ray observations 
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of nearby and bright AGN unveiled the co-existence of multiple 
absorbers/reflecting mirrors in the central regions, suggesting that 
the absorbers could be located at different scales (from within few 
tens of gravitational radii from the central nucleus t o outside the 
pc-scale torus) and coul d be in part inhomogeneous jRisaliti et al.l 
I2OO2I : lElvis et al. ||2004|) . Our view of the inner structure of AGN 
has been also modified by the recent finding of hard excesses 
in type 1 AGN, which can be modelled with the presence of a 
Compton-thick gas in the line-of-sig ht ( PDS 456, Reeves et al.l 
l2009t IH 0419-577, lTumer et al.ll20o"^ and lTatum et al. Il2012h . 



The variability of the column density of the X-ray absorbing 
gas (A^h) observed in a large number of AGN revealed that a 
significant fraction of the absorbing medium must be clumpy. 
The time-scales of these TVh variations, which can be directly 
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linked to the size and distance of the absorbing clouds, also 
provided valuable constraints on the size and location of this 
obscuring rnateria l from the central accreting black hole (see 
lRisalitietal]|2002l) . Rapid A^H variations have been discovered on 
time-scales from a few days down to a few hours for a limited but 
still increasing number of obscured or type 2 AGN : NGC 1365 



1980l): Gmgfl jAwaki etalJ Il99lh; BepnoSAX jRisalitil |2002| : 
Dadinall2007h : ASCA jTumer et al.|[T997l; IComastri et alill998ll : 



XMM-Newton\ Chandra 



and Rossi X-Ray 



Risaliti et al 



Elvis et al. 



^Matt et alj 

Timing Explorer (RXTE; [Rivers et al] 1201 ih . NGC 4507 is one 
of the brightest Seyfert 2s detected above 10 keV with the both 
BAT detector on board Swif t and INTEGRAL It is part of the 
58 months BAT catalogu fQ jTueller et al.ll2010l ; Baumgartner et 
2004), NGC 4151 IPu ccetti et al "2007 |), NGC 7582 al. 2012 ApJS submitted) and of the INTEGR AL AGN catalogue 



20051, l200l[2009l; 



Maiolino etallbOlOl), NGC 4388 



■Piconcelli et al. 2007: Bianchi et al.ii,2009 : Turner et alj|2000h and 
UCC 4203 ( iRisaliti et aljhoiOl) . Some of these extreme variations 
are effectively occultation events where the column density of 
the absorber changes from Compton-thin (A'h < 10^* cm~^) to 
Compton-thick (A'h > 10^'* cm^^). These A'^h variations unveiled 
that a significant fraction of such absorbing clouds must be located 
very close to the nuclear X-ray source and, more specifically, 
within the broad line region (BLR). 

However this picture is proven only for those few objects, 
which show extreme and rapid Nh variations. On longer time- 
scales (from months to years) A'h variability is a comm on property 
in local bright Seyfert 2 galaxies jRisaliti et al.ll2002l) . Due to the 
complexity of the measurements, when a change in the spectral 
shape is found it is hard to distinguish between A'h and photon 
index variations. The main open question for some of the detected 
variations is whether the spectral changes are indeed due to a 
variable circum-nuclear absorber or are due to variability of 
the intrinsic emission. In order to remove this degeneracy, high 
sensitivity and wide spectral coverage (to determine the contin- 
uum component) are needed. Variability of the X-ray absorbers 
is a common property of AGN. Indeed, in type 1 AGN (see 
[Turner & Miiieill2009l ; [Bianchi et al.ll2012h most of the observed 
spectral variability can be described with changes in the covering 
factors and ionisa tion states of the inner absorbers (N GC 4051, 
Miller etal.l bold: MCG-06-30-15, [Miller et al."2008': Mrk766, 



20-100 kcV I 



2 X 10" 



_er^ 



Miller et al 



l2007l : [Turner et all 20071 ; PDS 456, Behar et al. 2010). 



Furthermore, even if rare, o ccultation events ha ve been detected 
in type 1 AGN: M rk 766 jRisaliti et al J [201 ih and NGC 3516 
dTumer et all2008h . supporting the overall picture. 

The hypothesis of a clumpy structure for the absorbing 
"torus" has been recen t ly introduced in several theoretica l models 
jNenkova et al.[ [2002[ , [20081 : [Elitzur & Shlosman[ [20061 : [ElitzuJ 
[20121 and reference therein), where the torus consists of several 
distinct clouds, distributed in a soft-edge toru s. These models 
were originally base d on infrared observations jjaffe et aljliooi : 
[PonceletetaT] [20061) . showing an apparent si milarity between 
the IR emission of type 1 and type 2 AGN jLutz et alj [2004 
[Horst etalj|2006l) , but are now strongly supported by the short- 
term changes of the A'h of the X-ray absorbers. As discussed by 
[Elitzuj l l2012h . for a "soft-edged" toroidal distribution of clouds, 
the classification of type 1 and type 2 does not depend solely 
on the viewing angle: although the probability of a "unobscured 
view" of the AGN decreases when the line-of-sight is far from the 
axis, it is non zero. Furthermore, this model naturally accounts for 
A'h variability and in particular for occultations events due to the 
transition of a single cloud. 

NGC 4507 is one of the X-ray brightest (-F(2-io kcV) ~ 
0.6 - 1.3 X 10"^^ erg cm"^ s"^), and nearby (z = 0.0118) 
Seyfert 2 galaxies, with an estimated un absorbed luminosity of 
3.7 X 10*^ erg s~^ jComastri et alj[T998l) . It has bee n observed 
with all the major X-ray observatories: Einstein jKriss et al] 



_ Beckmann et alj2009l: 

[Malizia et all [20091 : [Bassani et all [2006h . NGC 4507 was also 
detected in the soft gamma-ray band with the OSSE experiment o n 
board the Compton Gamma Ray Observatory jBassani et al.ll995h . 

All of these X-ray observations revealed a hard X-ray spec- 
trum typical of a Compton-thin Seyfert 2: an X-ray continuum 
characterised by heavy obscuration and a strong Fe Ka line at 
6.4 ke V. The average measured column density is A^h ~ 6x 10^'^ 
cm "2 ( lRisaliti[|2002l) . A reflection component w i th a reflection 
fraction ranging from 0.7 to 2.0 jRisalitil [20021 : [Padind [2007h 
was measured with the BeppoSAX observations, while the high 
energy cut-off could not be constrained. A similar result has been 
obtained with the IN TEGRAL IBIS/ISGRI data (R = O.GlJ g; 
[Beckmann et alj2009l) . 

The soft X-ray spectrum is also typical of a Compton-thin 
Seyfert 2, with several emission lines from 0.6-3 keV range, due 
to ionised elements from O to Si, which require the presence 
of at least t wo ph otoionised media or a single stratified medium 
( [Matt et alj [2004h . A comparison between the observations 
performed with BeppoSAX and ASCA showed long-term A'h 
variability, which changes by a factor of 2, and also some possible 
variab ility of the intrinsic continuum ( [Risaliti et alj|2002 : RisalitH 
[2OO2I) . Altogether the emerging picture for NGC 4507 is of a 
complex and highly variable absorber, as seen in other bright 
Compton-thin Seyfert 2s. At least two absorbing systems are 
present: a Compton-thick reprocessor, responsible for the Fe Ka 
line at 6.4 keV plus the strong Co mpton reflected com ponent 
detected with BeppoSAX and RXTE ( [Rivers et al.1 [201 ih , and a 
variable Compton-thin absorber. A Chandra observation also 
provided a detection (at 99% significance) of a Fe XXV absorption 
line (at ~ 6.7 keV), whic h suggested the presence of an ionised 
absorber ( [Matt et al .l2004h . 

Here, we present the results of a Suzaku observation (of 
net exposure ~ 90 ks) of NGC 4507 and a comparison with 
previous X-ray observations of this AGN, which shows that below 
10 keV NGC4507 alternates from being in a transmission to a 
"reflection" dominated state. The Suzaku observation has been 
already presented in the statistical analy s is of 8 8 Seyfert 2 galaxies 
observed by Suzaku ( [Fukazawa et al.l [201 ih . investigating the 
properties of the Fe K comple x versus the amount o f absorption 
and luminosity of the sources. [Fukazawa et al.[|201 1[ report both 
the presence of a high column density absorber as well as a strong 
Fe K line complex {EW = 600 ± 30 eV). Here we present a more 
detailed analysis of the same observation, where we investigate the 
Fe K line complex allowing not only the Fe Ka but also the Fe K/3 
line parameters {E, a and /) to vary as well as a comparison with 
the previous X-ray observations. Finally, we also tested a new 
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model for the toroidal reprocessorFI jMurphv & Yaaoobll2009ll . 

The paper is structured as follows. The observation and data 
reduction are summarised in § 2. In § 3 we present the modelling 
of the broad-band spectrum obtained with Siizakii, aimed to assess 
the nature of the X-ray absorber, the amount of reflection and 
the iron K emission line. In § 4 we then compare the Suzaku 
with the previous XMM-Newton and BeppoSAX observations. 
Discussion and conclusions follow in § 5. Throughout this paper, a 
concordance c osmology with Hp = 71 km s^^ Mpc^^, f2A=0.73, 
and ^^=0.27 dSpergel et ai]|2003h is adopted. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Suzaku 

A Suzaku jMitsuda et al] l2007h observation of NGC 4507 was 
performed on 20th December 2007 for a total exposure time of 
about 103 ksec (over a total duration of ~ 178 ksec); a summary 
of the observations is shown in Table[T] Suzaku carries on board 
four co-aligned telescopes each w ith an X-ray CCD camera 
(X-ray Imaging Spectrometer; XIS iKovama et al.l 120071) at the 
focal pla ne, and a non imaging hard X-ray detector (HXD-PIN; 
iTakahash i et al. 1 12007|) . Three XIS (XISO, XIS2 and XIS3) are 
front illuminated (FI), while the XIS 1 is back illuminated; the latter 
has an enhanced response in the soft X-ray band but lower effective 
area at 6 keV than the XIS-FI. At the time of this observation 
only two of the XIS-FI were still operating namely the XISO 
and XIS3. All together the XIS and the HXD PIN instruments 
cover the 0.5-10 keV and 14-70 keV bands respectively. The 
cleaned XIS event files obtained from version 2 Suzaku pipeline 
processing were processed using HEASOFT (version v6.6.3) and 
the Suzaku reduction and analysis packages applying the standard 
screening for the passage through the South Atlantic Anomaly 
(SAA), elevation angles and cut-off rigiditjQ. The XIS data were 
selected in 3 x 3 and 5x5 editmodes using only good events with 
grades 0, 2, 3, 4, 6 and filtering the hot and flickering pixels with 
the script sisclean. The net exposure times are 87.7 ksec for each 
of the XIS and 92.9 ksec for the HXD-PIN. 

The XIS source spectra were extracted from a circular region 
of 2.6' radius centered on the source, while background spectra 
were extracted from two circular regions of 2.6' radius offset 
from the source and the Fe 55 calibration sources, which are in 
two corners of CCDs. The XIS response (rmfs) and ancillary 
response (arfs) files were produced, using the latest calibration 
files available, with the ftools tasks xisrmfgen and xissimarfgen 
respectively. The spectra from the two FI CDDs (XIS and XIS 
3) were combined in a single source spectrum (hereafter XIS-FI) 
after checking for consistency, while the BI (the XISl) spectrum 
was kept separate and fitted simultaneously. The net 0.6-10 keV 
count rates are: (0.128 ± 0.001) counts s"\ (0.131 ± 0.001) 

^ http://www.mytorus.com/ 

3 The XIS 2 failed on November 2006 

^ The screening filters all events within the SAA as well as with an Earth el- 
evation angle (ELV) < 5° and Earth day-time elevation angles (DYE_ELV) 
less than 20°. Furthermore, we excluded also data within 256s of the SAA 
from the XIS and within 500s of the SAA for the HXD. Cut-off rigidity 
(COR) criteria of > 8 GV for the HXD data and > 6 GV for the XIS were 
used. 



counts s-\ (0.128 ± 0.001) counts s~^ for the XIS 0, XIS3 
and XISl respectively. Data were included from 0.6-10 keV for 
the XIS-FI and 0.6-8.5 keV for the XISl chip; the difference on 
the upper boundary for the XISl spectra is because this CCD is 
optimised for the soft X-ray band. We also excluded the data in 
the 1.6-1.9 keV energy range due to calibration uncertainties. The 
XIS FI (BI) source spectra were binned to 1024 channels and then 
to a minimum of 100 (50) counts per bin, and statistics have 
been used. 

The HXD-PIN data were reduced following the latest Suzaku 
data reduction guide (the ABC guide Version 4.(0), and using 
the rev2 data, which include all 4 cluster units. The HXD-PIN 
instrument team provides the background event file (known as the 
"tuned" background), whi ch accounts for the i nstrumental Non 
X-ray Background (NXB; iKokubun et al. ||2007|) . The systematic 
uncertainty of this "tuned" background model is believed to be 
±1.3% (at the io level for a net 20 ks exposure). We extracted 
the source and background spectra using the same common good 
time interval, and corrected the source spectrum for the detector 
dead time. The net exposure time after screening was 92.9 ks. We 
simulated a spectru m for the co smic X-ray background counts 
adopting the form of lBoldJ ( Il987h and lOruber et al.lll999l) and the 
response matrix for the diffuse emission; the resulting spectrum 
was then added to the instrumental one. 

NGC 4507 is detected up to 70 keV at a level of 23.6% 
above the background corresponding to a signal-to-noise ra- 
tio S/N ~ 50. The net count rate in the 14-70 keV band is 
0.122 ± 0.002 counts s"^ (~ 11300 net counts). For the spectral 
analysis the source spectrum of NGC 4507 was rebinned in order 
to have a signal-to-noise ratio of 10 in each energy bin. A first 
estimate of the 14-70 keV flux was derived assuming a single 
absorbed power-law (T = 1.8 ± 0.3) model. The 14-70 keV flux 
is about ~ 9.2 X 10~^^erg cm^'^ s^^ and the extrapolated flux 
in the Swift band (14-195 keV) is ~ 1.8 x 10""' erg cm"^ s"\ 
comparable to the flux reported in the BAT 58 months catalog 
(~ 1.9 X lO'^^erg cm~^ s"\ Baumgartner et al. 2012 ApIS 
submitted). 



2.2 XMM-Newton 

XMM-Newton observed NGC 4507 on the 4th January 2001 
(see Table[T) with a total expos ure time o f about 46.2 ksec. This 
observation was presented by iMatt et ah I l l2004 . Since we are 
mainly interested in a comparison with the Suzaku observation and 
since a detailed analysis has been already published, we focused 
only on the EPIC-pn data. The EPIC-pn camera had the medium 
filter applied and was operating in Full Frame Window mode. The 
XMM-Newton data were processed and cleaned using the Science 
Analysis Software (SAS version 10.0.2) and analysed using stan- 
dard software packages and the most recent calibrations. In order 
to define the threshold to filter for high-background time intervals, 
we extracted the 10-12 keV light curves and filtered out the data 
when the light curve is 2cr above its mean. This screening yields 
net exposure time (which also includes the dead-time correction) 
of ~ 32 ksec. The EPIC pn source spectrum was extracted using 
a circular region of 35" and background data were extracted using 

^ http://heasarc.gsfc.nasa.gov/docs/suzaku/anaIysis/abc/ 
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two circular regions with a radius of 30" each. Response matrices 
and ancillary response files at the source position were created 
using the SAS tasks arfgen and rmfgen. The source spectrum was 
then binned to have at least 50 counts in each energy bin. 



2.3 The Swi/if-BAT 

The Swift-BKi spectrum and the latest calibration response file (the 
diagonal matrix: diagonal.rsp) were obtained from the 58-month 
survey archive; the data r eduction procedure o f the eight-channel 
spectrum is described in iTueller et alj feOld) and Baumgartner 
et al. (2012 ApJS submitted). The net count rate in the 14-195 
keV band is (2.69 ± 0.05) x 10~^ counts s~^ (corresponding to 
Fi4-i95kcv ~ 1.9 X 10-i°ergcm-2 s"'). 



3 SPECTRAL ANALYSIS 

All the models we re fit to the data using standard software packages 
(XSPECver. 12.7.0 lArnaudll996l) and including Galactic ab sorption 
(A^H.Gai = 7.23 X 10^°cm~^: lbickev & Lockmanlll990h . In the 
following, unless otherwise stated, fit parameters are quoted in the 
rest frame of the source and errors are at the 90% confidence level 
for one interesting parameter (Ax^ = 2.71). 



3.1 The 0.6-150 keV continuum 

For the analysis we fitted simultaneously the Suzaku spectra 
from the XIS-FI (0.6-10 keV), the XIS 1(0.6-8.5 keV) HXD-PIN 
(14-70 keV) and the Swift-hKi spectrum (14-150 keV). We set 
the cross-normalization factor between the HXD and the XIS-FI 
spectra to 1.16, as recommended for XIS nominal observation 
processed after 2008 Jul30 (Manabu et al. 2007; Maeda et al. 
2008), while we left the cross-normalisation with the Swi/r-BAT 
spectrum free to vary. 

We initia lly tested the best-fit continuum model 
presented in iMatt et all j2004h for the XMM-Newton 
data, which was of the mathematical form: F{E) = 
wabs X (zwabs x powl + pexrav + pow2), where POWl 
is the absorbed power-law, PEX RAV is the XSPEC mod e l for a 
Compton-reflected component jMagdziarz & Zdziarskil Il995h . 
POW2 is the soft scattered power-law continuum, which is ab- 
sorbed only by the local Galactic absorber (wabs) and zwabs is the 
intrinsic absorber. For the PEXRAV component we allowed only 
its normalisation to vary, while we fixed the high energy cutoff 
to 200 keV, the amount of reflection R = Q/2n — 1 and the 
inclination angle i to 30° , as adopted by the previous work on the 
XMM-Newton data. 

As shown in Fig. [T] the residuals to our baseline continuum 
model at the energy of the Fe K band clearly reveal the presence of 
a strong narrow core at the expected energy of the Fe Kq (6.4 keV), 
as well as a strong Fe K/3 (~ 7.06 keV). We then included two 
narrow Gaussian lines to account for Fe Kq and Kfi emission lines 
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Figure 1. Data/model ratio between the XIS data (XIS-FI, black filled 
squares in the electronic version; XIS-BI red open triangles in the elec- 
tronic version) and the continuum model, showing the iron line profile. The 
two vertical dashed lines correspond to the rest-frame energies of the Fe Ka 
and K/3 emission lines at 6.4 keV and 7.06 keV respectively. 



(see §3.2). Initially we fixed the energy of the Fe K/3 emission line 
to 7.06 keV, we tied its width to the width of the Fe Ka emission 
line and its normalisation to be 13.5% of the Fe Ka emission line 
jPalmeri eTal. .2003) . The inclusion of the lines improves the fit 
by Ax^ = 1511 for 3 degrees of freedom (x^/d.o.f. =723.7/384); 
statistically the fit is still unacceptable, with most of the remaining 
residuals being at i5 < 2 keV. Following the results of previous 
X-ray studies of NGC 4507, which showed a remarkably steep 
X-ray emission below 2 keV ac companied by several soft X-ray 
emission lines ( lMattetal.ll2004h . we then allowed the scattered 
component (P0W2) to have a different photon index (F; see Table 
|4]l with respect to the primary power law and we added several 
narrow Gaussian emission lines, which even at the Suzaku CCD 
resolution are clearly visible (see Fig.|2j. 

We detected 6 strong soft X-ray emission lines (from Ne, Mg 
and S; see Table[2j and although relatively simple, this model al- 
ready provides a good fit (x^/d.o.f. = 416.0/372). However, we 
note that the soft power-law component is unusually steep (F — 
3.8 ± 0.2). The steep photon index could be due to the presence of 
other weak emission lines, which are unresolved at the XIS-CCD 
resolution and could be related either to the photoionised emitters 
responsible for the str ongest emission lines as seen in Compton- 
thin Seyfert galaxies jBianchi et al.l |2006| ; iGuainazzi & Bianchil 
l2007h or to the presence of additional emission from a collisionally- 
ionised diffuse gas. Although we investigated different physical 
scenarios (see ^3.2\ , we note that the different models tested for 
the soft X-ray emission did not strongly affect the results of the 
hard X-ray energy band. 



3.2 The soft X-ray emission 



Although our primary aim is the analysis of the hard X-ray 
emission, we investigated both photoionised and collisionally 
ionised plasmas as sources for the soft X-ray emission lines; to 
this end we fitted the 0.6-150 keV spectra replacing in turn the 
http://wvvw.astro.isas.jaxa.jp/suzaku/doc/suzakumemo/suzakumemo-2008-06.pdfl Gaussian emission lines with either an additional thermal com- 



° http://www.astro.isas.jaxa.jp/suzaku/doc/suzakumemo/suzakumemo- 
2007-1 1. pdf; 
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Table 1. Summary of the observations used: Observatory, Epoch, Instrument, and Net exposure times. The net exposure times are after the screening of the 
cleaned event files. 



Mission 


DATE 


Instrument 


T(nct) (ks) 


Suzaku 


2007-12-20 


XIS 


87.7 


Suzaku 


2007-12-20 


HXD-PIN 


92.7 


XMM-Newton 


2001-01-04 


EPIC-pn 


32.3 


BeppoSAX 1 


1997-12-26 


MECS 


49.79 


BeppoSAX 1 


1997-12-26 


PDS 


26.97 


BeppoSAX 2 


1998-07-02 


MECS 


31.41 


BeppoSAX 2 


1998-07-02 


PDS 


16.95 


BeppoSAX 3 


1999-01-13 


MECS 


41.33 


BeppoSAX 3 


1999-01-13 


PDS 


20.08 



Table 2. Summary of the strongest soft-X-ray emission lines as detected in the Suzaku spectra. The energies of the lines are quoted in the rest frame. Fluxes 
and possible identifications are reported in column 2 and 3. The observed EW are reported in column 4 and they are calculated against the total observed 
continuum at their respective energies. In column 5 the improvement of fit is shown with respect to the continuum model; the value for the model with no soft 
X-ray lines is x^/d.o.f. = 723.7/384. Finally in column 6 we report the lab energy for the detected fines. 



Energy 
(keV) 
(1) 


Flux 

(IQ-f'phcm-^ s-i) 
(2) 


ID 

(3) 


EW 
(eV) 
(4) 


Ax" 
(5) 


ELab 

(keV) 
(6) 


90+°-0^ 
"■^^-0.01 


2&.6tll 


Ne IX He-a 




164.2 


0.905(f); 0.915(i);0.922 (r) 






Ne X Lya 




21.5 


1.022 


, 9, +0.01 
'■'^^ -0.01 




Ne X Ly/3 


471^3 


20.8 


1.211 (r) 


^■■'°-o.oi 


■'■■^-i.i 


Mg XI Heo 


'"-15 


69.4 


1.331(f); 1.343(i); 1.352 (r) 


2 41+0.03 
^■^'-0.05 


2 0+0.9 


S XIV Ka 




12.5 


2.411 


, 7r,+0.03 
-"■'"-0.03 


^ 1+0.8 

^■'-0.8 


CaKa 


co+20 
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ponent (MEKAL model in XSPEC. lMewe et aljri985l) or a grid of 
photo ionised emission model generated by XSTAR Jkallman et al.l 
l2004h . whichi assumes a F ~ 2 illuminating continuum and a 
turbulence velocity of — 100 km/s. 

We found thiat neithier a single photoionised emission model 
or a thermal component provide an acceptable the fit (x^/d.o.f. 
=592.3/382 and x^/d.o.f. = 584.7/382 for the XSTAR and MEKAL 
component with respect to the model with no soft X-ray lines 
X^/d.o.f. =723.7/384, respectively), and strong residuals are 
present below 2 keV; furthermore, for both these models we found 
a steep F for the the soft power-law component (F = 3.2 ± 0.2). 
Thus we tested for the soft X-ray emission a composite model 
consisting of: a collisionally ionised emitter, a photoionised 
plasma and a soft power-law component. We found that this 
model is now a better representation of the observed emission 
(X'^/d.o.f. =510. 1/380). The photoionised emitter has an ionisation 
parametefl of log^ = LOG^o Jl erg cm s~^, while the thermal 
component has a temperature of kT — 0.75 1 g 1)5 ^'^'^ 
model gives a total observed (i.e. corrected only for the Galactic 
absorption) 0.5-2 keV flux of ~ 4.8 x lO^^^erg cm~^ s"^ in 
this energy range the relative contribution of the collisionally and 
photoionised emitters are Fg ^Ifi ~ 9.7 x 10""erg cm"^ s~^ 



^ The ionisation parameter is defined as ^ = Lion/ncR^, where L is the 
ionising luminosity from 1-1000 Rydberg, R is the distance of the gas, ne 
is its electron density. 



and Fo fS^ ~ 1.3 x lO^^^erg cm"^ s"^ respectively. We note 
that the limited spectral resolution of the CCD spectra prevents 
us from deriving definitive conclusion on the relative importance 
of these two emission components, furthermore the photon index 
is still relatively steep F — 3.0 ± 0.3, suggesting that there is 
still a possible contribution from unresolved emission lines or a 
collisionally ionised emission (e. g. a thermal component). 
Recently we obtained a deep (~ 130 ks) XMM-Newton-RGS 
observation of NGC 4507, which provided the best soft X-ray 
spectrum so far for NGC 4507. The properties of the soft X-ray 
emission are discussed in more detail in a companion paper 
describing the XMM-A^ew ?on-EPIC and XMM-A'ewton-RGS data 
jMarinucci et al. 1 1201 2bl . Wang et al. private communication). 
Briefly the RGS data unveiled that the soft X-ray emission is 
indeed dominated by emission lines, and that cannot be explained 
with a single photoionised or thermal component. We stress again 
that the different models tested for the soft X-ray emission did not 
strongly affect the results of the hard X-ray emission, which is our 
primary focus in this paper. 



3.3 The Fe K emission line complex and the high energy 
spectrum 

We then considered the hard X-ray emission of NGC 4507 using 
for the soft X-ray emission the simple phenomenological model of 
a scattered power-law component and 6 Gaussian emission lines. 
The spectrum was then parametrized with a model of the form: 
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Figure 2. Data/model ratio between the XIS data (XIS-FI, black filled 
squares in the electronic version; XIS-BI red open triangles in the electronic 
version) and the continuum model (see § 3.1), showing the soft X-ray lines. 



F{E) — wabs X (zwabs x powl + pexrav + Fe Ka + Fe K/3 + 
pow2 + 6 GAom), where the ratio of Fe Ka and Fe K/? intensities 
was initially iixed at 13.5% and GAem are the soft X-ray emission 
lines. As previously described, the photon index of the scattered 
power-law component (P0W2) was left free to vary independently 
from the primary power-law component (POWl). This model pro- 
vides a good description of the continuum (x^/d.o.f. = 416.0/372); 
an intrinsic column density of A^h = (8.4 ± 0.5) x 10^'' cm^^ 
is required, the photon index of the primary absorbed power 
law is r = 1.83 ± 0.04 and the intensity of neutral reflection 
component is Apexrav = (1.2 ±0.1) x 10^^ photons cm^^s^^ 
(corresponding to _R ~ 1.6). 

Taking into account the high statistics of the present data and 
the strength of both the Fe K lines, we then left free to vary the 
centroid energy and intensity of the Fe K/3 line and we found a 
statistically similar best-tit (x^/d.o.f =409.5/370). For the Fe Ka 
line core we obtained E = QAQ^tol^ keV, a = 35 ± 10 eV 
(corresponding to a FWHM = 3860 ± 1100 km s"^) and 
EW — 490 ± 30 eV (with respect to the observed continuum). 
For the corresponding Fe K/3 we obtained a centroid energy of 
E = 7.07 ± 0.02 keV and /fck/3 = 9.50ti:™ x lO^'^photons 
cm~^s~^, which corresponds to a /fcKq/ZfcK^ ratio to about 
~ 17%. We also checked the accuracy of the energy centroids 
and line widths using the ''''Fe calibration sources located on two 
corners of each of the XIS chips, which produce lines from Mn Ka 
(Kal and Kq2 at 5.899 keV and 5.888 keV respectively). From 
measuring the lines in the calibration source, we find no major 
energy shift or residual broadening (E — 5.903 ± 0.003 keV, 
a < 15 eV). Thus the apparent broadening of the Fe Ka emission 
line is not due to calibration uncertainties, however, upon the 
inclusion of a possible Compton shoulder the Fe Ka is unresolved 
(see below and §4.1). 

The parameters der ived for the Fe K line c omplex are similar 
to the values reported in lFukazawa et al. l|); we note however 
that in our modelling the energy centroids and normalisations 
of all the lines are left free to vary. In the upper panel of Fig. [3] 
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Figure 3. Upper panel: 68%, 90%, and 99% confidence contours of en- 
ergy centroid of the Fe Ka versus the Fe K/3 emission fines (rest frame). 
Lower panel; 68%, 90%, and 99% confidence contours of the intensity of 
the Fe Kq versus the Fe K/3 emission line; the dashed fine corresponds to 
the Fe K/3/Fe Ka line ratio of 0.135. 



we show the 68%, 90%, and 99% confidence contours of the 
narrow Fe K/3 centroid energy versus the centroid energy of the 
Fe Ka, while in the lower panel we show the coiTesponding 
contours for the line intensities. We note that contamination to the 
Fe K/? from a possible Fe XXVI {E = 6.97 keV) is negligible, 
indeed as can be seen in Fig. [3] (upper panel) the contours of 
the energy centroid of the Fe K/3 are fairly symmetric and not 
elongated toward lower energies. Although the ratio of Fe K/3 
and Fe Kq intensities is consistent within the errors (at 99%, see 
Fig. [3] lower panel) with the 13.5 % value, as expected for low 
ionisation Fe jPalmeri et al.l |2003|) . it is marginally higher than 
the theoretical value for low ionisation Fe; such a high value 
of the Fe K/3/Fe Ka could be indicative t hat the Fe ionisation 
state could be as high as Fe IX. In particular IPalmeri et alj ( |2003|) 
showed that while for Fe I both theoretical and experimental 
values of this ratio lie in the 12-13.5% range, for higher ionisation 
states this ratio can be higher and for Fe IX it can be as high as 17%. 

We note that the parameters of the continuum (see Table |4ll 
and in particular of the reflection component are all well con- 
strained. In Fig. 151 we show the confidence contours between 
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the normalisation of the reflection component and the intrinsic 
photon index (F), obtained allowing the cross normalisation 
factor of the HXD over the XIS-FI spectrum to vary. The best 
fit value of the intensity of the neutral reflection component is 
A'^PBXRAv = 1.2 ± 0.1 X 10"^ photons cm"^ s"\ which 
corresponds to i? = 1.9 ± 0.2 ; we note that this component 
dominates the spectrum below ~ 7 keV (see Fig.[5j. 

We found also evidence for the Fe Ka Compton shoulder 
(E = 6.23 ± 0.08 keV), which is significant at the 99.8% 
confidence level, according to the F-test (Ax^ = 11.5 for 2 dof, 
X^/d.o.f. =398.0/368). The ratio between the Compton shoulder 
intensity and that of the Fe Ka is ~ 9%, which together with the 
strong Compton reflected c omponent confirms the presence of a 
Compton-thick reprocessor ( lMatlll2002l : lYagoob & Murphvll201lh . 

We also tested for the presence of emission lines from Fe XXV 
(E = 6.7 keV) and Fe XXVI {E = 6.97 keV) adding to the model 
two narrow Gaussian emission lines; the former is detected at 
E ~ 6.73 keV (Ax^ = 8.1, xVd-o.f. =389.0/366; see Tabled 
with an EW~ 26 eV, while for the Fe XXVI emission line we can 
place an upper limit on its flux to Jpg xxvi ^ ^-^ ^ 10^'' photons 
cm^^s^^ (or EW < 26 eV). Finally, the inspection of the 
residuals left by this model unveils the presence of an emission 
line like feature at ~ 7.5 keV, suggesting the presence of emission 
from Ni Ka. Thus we included an additional narrow Gaussian 
line and we found that the parameters of this additional line are 
consistent with Ni Ka (E = 7.50 ± 0.08 keV, Ax^ = 6.7 for 
2 d.o.f.). This model now provides a good description of the 
broadband X-ray spectrum of NGC 4507 (xVd-o.f. =383.2/364) 
and no strong residuals are present (see Fig.[5](. After the inclusion 
of these additional narrow Gaussian lines the Fe Ka line is now 
unresolved with cr < 30 eV corresponding to FWHM < 3300 
km s~^. The Fe Ka width measured with Suzaku is in agreement 
with the Chandra HETG measurement of cr < 26 eV iMatt et a ll 
l2004h . and suggests an origin in the torus. We note that there are 
no residuals left in the Fe K region for a possible strong underlying 
broad component {EW < 120 eV); furthermore, leaving the width 
of the other emission lines free to vary does not improve the fit 
(Ax' = 1). 

In order to understand if the high value of the reflection 
component could be due to the adopted model for the reflec- 
tion component, we tested the COMPPS model developed by 
iPoutanen & SvenssorJ ( Il996l) . This model includes the processes 
of thermal Comptonisation of the reflected component, which is not 
included in the PEXRAV model. We tested both a slab and a spher- 
ical geometry for the reflector. We found that both these models 
provide a good description of the observed emission (x' /d.o.f. = 
400.1/361 and xVd-o.f. = 394.5/360 for the slab and spherical 
geometry respectively) with no clear residuals. In both these sce- 
narios the reflection fraction is similar to the one measured with 
the PEXRAV model {R > 1.6 in both cases). 

The amount of reflection is cons i stent with the Beppo SAX 
measurements of NGC 4507 l lRisalitil |2002| ; iDadinal |2007|) . for 
whi ch the authors report a reflection fraction ranging from 0.7 to 
2.0 l lRisaUti|[2003) . while it is remarkably higher th an the value 
repor ted from a RXTE measurement (R — 0.4 ± O.l. lRivers et al.l 
I2OIII) . The apparent discrepancy between these measurements 
could be ascribed to the combination of different effects; among 
them, variability of the primary continuum and of the amount of 
absorption (see §4). In particular, the Suzaku measurement appears 




Figure 4. Contour plot showing the 68%, 90%, and 99% confidence levels 
of the normalisation of the reflected component versus F. Note that the 
contours have been produced allowing all the other parameters of the best- 
fit model to vary, including the cross-normalisation between the Suzaku XIS 
and HXD. 



to be, at a first glance, consistent with the scenario proposed for 
the BeppoSAX observations, where the increase of the reflection 
fraction was ascribed to a Compton-reflected component remaining 
constant despite a drop in the primary power-law flux. We note 
that the RXTE observations were performed in two campaigns one 
in 1996 and one in 2003, with 94% of the total good exposure 
time being from the 1996 campaign. During these observations 
the observed 2-10 keV was a factor of two higher than during 
the Suzaku one, and thus the lower reflection fraction could be in 
agreement with this scenario. 



However, as already suggested by iRivers et"aL I feOIll) . the 
scenario could be far more complex and also indicative that the 
reflection component normalisation is responding to a different 
past illuminating flux. We note however that not only do these 
works assume different values for the inclination angle, which 
could affect the measurement of the reflection fraction, but that 
also the model itself which is adopted for the reflected component 
is not flawless, indeed it assumes that the reflector is a semi-infinite 
slab and also its density is assumed to be infinite. More importantly 
the energy band and spectral resolution of the observations have a 
strong impact on the measurement of the continuum parameters; 
for example, the RXTE observations have a lower resolution at the 
energy of the iron line with respect to the BeppoSAX and Suzaku 
ones, and this has a strong impact on the measurement of the Fe K 
line/edge properties as well as on the measurement of the amount 
of reflection and the intrinsic F. The amount of reflection also 
strongly depends on the model assumed for the X-ray absorber, 
indeed including the effect of the Compton-down scattering would 
increase the normalisation of the primary emission and not of 
the reflected component and thus lower the value of the reflected 
fraction. 

A more detailed description of the variability of this source is 
presented in §§4 and 4.1, where we compare the historical X-ray 
spectra obtained for this source, showing that both the amount of 
absorption and possibly the primary continuum are varying, while 
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Figure 5. Suzaku and Swift 0.6-150 keV data and best-fit model (in the 
electronic version: XIS-FI, black; XISl, red; HXD-PIN green; Swift-B AT 
blue) of NGC 4507; data have been rebinned for plotting purposes. The 
upper panel shows the data and model (F ~ 1.8; Njj ~ 8 X 10^^ cm~^see 
Table|4]and §3.3), fitted over the 0.6-150 keV band. The lower panel shows 
the data/model ratio to this model. 



the reflected component and the A'h of a distant reprocessor remain 
rather constant. We present a summary of the main parameters of 
the X-ray emission that can be derived assuming both the standard 
models and th e new model for a toroida l reprocessor, i.e. the MY- 
TORUS model l lMurphv & Yaqoobllioogh . 



4 EVIDENCE FOR A VARIABLE ABSORBER 

In Fig.|6]we compare the XMM-Newton (red) and the Suzaku XIS 
& HXD (black) data; a clear difference in the curvature is present 
between 4 and 8 keV, which is most likely due to a change in the 
amount of absorption of the primary radiation. 

To test this hypothesis we applied the Suzaku best-fit model 
to the XMM-Newton spectrum, allowing the Fe Ka emission line 
parameters as well as all the continuum parameters free to vary 
(see Table |4). We found a statistically acceptable fit (x'^/d.o.f. = 
633.4/503), which unveiled that the main difference between these 
two observations can be explained with a lower column density 

5.0 ± 0.3 X 10^^ cm-^and 
8.2 ± 0.6 X 10^^ cm-2). The normahsation of 



of the intrinsic absorber {N^^^ 



the primary power-law component also varied between the two 
observations as well as the amount of reflection; however, we 
note that there could be some degeneracy between the slope of 
the primary power-law component and the amount of reflection 
when lacking a simultaneous high energy observation. We note 
that without allowing the column density to vary, we could not 
reproduce the different 2-10 keV spectral curvature observed 
during the XMM-A'ewto;i observation. 

To break this degeneracy and better understand the variability 
of NGC 4507 we reanalysed the 3 BeppoSAX observations of 
NGC 4507 (hereafter SAXl, SAX2 and SAX3 see Table 1) and, 
since we are mainly interested in the hard X-ray emission and 
variability of the spectral curvature, we considered only the MECS 
and PDS data in the 2-10 keV and 15-200 keV energy range 
respectively. We adopted the best-fit model of the Suzaku data and 
we fixed the components responsible for the soft X-ray emission 
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Figure 6. Comparison between the X-ray emission measured with Suzaku 
in 2007 (black data points in the electronic version) and XMM-Newton (red 
data points in the electronic version) in 2001. The continuum model is a 
standard model; composed of a primary power-law component transmitted 
trough a neutral absorber (labelled with T), a scattered power-law compo- 
nent (labelled with S) and a Compton reflected component (labelled with 
R and modelled with PEXRAV). For the simultaneous fitting we allowed to 
vary both the A^h and the normalisation of the primary power-law compo- 
nent. 



to the Suzaku values. We found that a simple change of the amount 
of X-ray absorption and the photon index cannot explain the 
observed variability. We then allowed also the nomalization of 
the primary continuum and the Fe Ko line parameters to vary 
while constraining the normalization of the reflected component 
to scale with the primary continuum (i.e. we fixed the ratio R, 
between the primary and reflected component to the one measured 
during the Suzaku observation). This model represents a situation 
where the reprocessor responsible for the reflected component 
responds to the variability of the primary continuum and thus it 
implies that this absorber should be close to the primary X-ray 
source. This model did not provide a good fit to SAXl, SAX2 
observation (xVd-o.f. =285.9/145, xVd-o.f. = 159.7/108), while 
it is a statistically acceptable fit for the last BeppoSAX observation 
(xVd-o.f. =110.9/100), during which NGC 4507 was in a state 
similar to the Suzaku observation. 

Upon allowing also the ratio between the intensity of the re- 
flection component and the primary continuum free to vary (i.e. 
the parameter R) the fits were acceptable and we found x'^/d.o.f. 
= 159.6/144, xVd-o.f. = 101.5/107 and xVd-o.f. = 108.6/99 for 
the SAXl, SAX2 and SAX3 observations, respectively. The pa- 
rameters of these best-fits are reported in Table|4j we no te that they 
are in a greement with the r esults previous ly presented in lMatt et al.l 
j2004l). lRisalitil ||2002|) and lOadinal (l2007h . We found that the Fe K 
line emission complex is rather constant and also note there is no 
evidence for variability of the intensity of the reflection component 
with the BeppoSAX observations. 



4.1 A more physical model 

This simple test, as described above, shows us that the variability 
properties of NGC 4507 are more complex than a simple variation 
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of the amount of absorption. We also note that the intrinsic 
photon index as well as the Fe K emission line complex are 
not variable, while variations are present in the column density 
and intensity of the continuum level and thus in the ratio of the 
reflection component versus the primary continuum. However 
the absolute flux of the reflection component is consistent with 
being constant. This in turn tells us that there is a rather stable 
reprocessor, which is responsible for the Fe K emission lines and 
the reflected component and which does not appear to respond to 
the variability of the primary continuum. This could be indicative 
of distant reprocessor, which does not respond to the variability of 
the primary continuum. Alternatively as we will discuss later this 
could indicate a clumpy absorber where the overall distribution 
of clouds remains rather constant. Given the limitations of the 
PEXRAV model, already outlined above (i.e. the geometry and 
density assumed for the reflector), this simple model does not 
allow us to derive strong constraints on the true nature of the 
absorber. Furthermore, by adopting this non physical model the 
temporal properties of the reprocessor (i.e. the variability of the 
amount of line-of-sight absorption and Compton reflection) could 
be highly uncertain and degenerate with respect to the variability 
of the primary continuum. Finally, the column densities of the 
X-ray absorber are in the range where the correction for the 
Compton-down scattering starts to be important. Thus we first 
included an additional absorber (CABS model in XSPEC) to 
account for the effect of the Compton-down scattering. We found a 
similar trend (as the one reported in Table 4) in the normalisations 
of the primary power-law components and thus in the intrinsic 
2-10 keV luminosities, albeit with a larger spread. 

Therefore we decided to reanalyse the available spec- 
tra using the most rece nt model for the toroidal reprocessor 
jMurphv & Yaqoobll2009l) . which correctly accounts for the emis- 
sion expected in transmission (hereafter zeroth order continuum) 
and reflection and also includes the expected Fe K emission lines 
(Fe Ka, Fe K/3 and the Compton shoulder). The calculations at 
the basis of this new model are all fully relativistic and valid 
for A'h in the Compton-thin and Compton-thick regimes. This 
new model assumes a uniform and essentially neutral toroidal 
reprocessor with an opening angle of 60° with respect to the axis 
of the system, while we note that the PEXRAV model assumes a 
disc/slab geometry for the reflector and thus the parameters derived 
from this model, such as the covering factor, can not be directly 
related to a covering factor of the putative torus as well as to the 
line-of-sight column density. In summary in the MYTORUS model 
all the different continuum components (reflected and transmitted 
components) and the fluorescent emission lines are all treated 
self-consistently and can thus be all directly related to the key 
parameters of the matter from which they originate. By adopting 
this model we were able to determine which component dominates 
in each energy band (reflected or transmitted components), assess 
their variability properties and thus better understand the global 
distribution of the absorber. 



4.1.1 A new implementation of the Mytorus model applied to the 
Suzaku observation 

The standard MYTORUS model, developed for XSPEC, is composed 
of three tables of reprocessed spectra calculated assuming that 
the input spectrum is a power law. These tables correspond to 
the main model components expected from the interaction of 



the primary power-law component with a reprocessor that has a 
toroidal geometry: the distortion to the zeroth-order (transmitted) 
continuum (MYtorusZ), the reflected continuum (MYtorusS), and 
the Fe Kq, Fe K/3 emission-line spectrum (MYtorusL). MYtorusZ 
is a multiplicative table that contains the pre-calculated transmis- 
sion factors that distort the incident continuum at all energies due 
photoelectric absorption. 

We first applied thi s toroidal-reprocessor model (MYTORUS 
iMurphv & YaqoobI |2009|) to the Suzaku observation. The model 
setup is the following: 

PHABSX (APEC -I- APEC + Agoft XZPOWERLW -I- 6 GAem 

-I- Fexxv -I- Ni Ka -I- MYtorusZ xzpowERLW -f Arx MY- 
torusS -I- Alx gsmooth X MYtorusL) 

We also included the Fe XXV and Ni Ka emission lines 
as well as a soft power-law component (^soft x zpowerlw) that 
represents scattering off optically-thin ionised gas (warm or hot), 
which are not included in the MYTORUS model. For the soft 
X-ray emission we kept the 6 Gaussian emission lines (GAcm; see 
Table 2) and we also included two thermal emission components, 
which allowed us to tie the photon index of the soft power-law 
component to the primary one as expected from scattering off 
optically-thin ionised gas. The normalisations of the MYTORUS 
components and of the soft power-law component are all tied 
together, while the value of the relative normalisations are included 
in the factors Asoft, and Al, which are in turn the relative 
normalisations of the soft power-law component, of the reflected 
component and of the emission lines. Since we expect the size 
scale of the scattering/reflecting and line emitting regions to be 
similar, we initially set Ab. ~ ~ 1 to be equal and we 
set them to 1. We note that we can not interpret any difference 
between these two factors as a difference in the size scales of 
these zones since these constants also include the effects of the 
transfer functions of the reflected continuum and line spectrum. 
The gsmooth component is the broadening of the Fe K emission 
lines and it is actually composed of two Gaussian convolution 
components, one is the actual broadening of the Fe Ka emission 
line while the second component accounts for the weak residual 
instrumental broadening as measured with the calibration sources 
{a < 15 eV, with a \/£ dependence). 

An inspection of the Suzaku spectra shows that above 8 keV 
we are dominated by the primary component, transmitted through 
the reprocessor, which is mainly constrained by the high-energy 
excess above 10 keV, while the reflected component dominates 
below 8 keV, where the strong emission lines from the Fe Ka, 
Ca Ka and S XIV Ka are present (see Fig. 7). This is analogous 
to what is observed adopting the old PEXRAV model (see Fig. [5] 
upper panel) and it is also in agreement with the variability of the 
normalisation of the primary power-law component as suggested 
in section §4. 

We note that in the standard configuration of the MYTORUS 
model, we cannot account self-consistently for the requirement of 
a strong transmitted component emerging at higher energies, the 
intense Fe K emission lines and a dominant reflected component 
below 8 keV. In particular, if we adopt the standard toroidal 
geometry with the inclination angles (between the axis of the 
torus and the observer's line of sight) of the two components 
tied together and if we do not allow the normalisations of the 
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Figure 7. The 0.5-70 keV Suiaku spectra when fitted with the MYTORUS 
model (in the electronic version: XIS-FI, black; XISl, red; HXD-PIN 
green). The dashed lines represent the reflected component (marked with 
R), the dotted lines are the transmitted component (marked with T) and the 
scattered component (marked with S) from ionised gas. 

reflected component to vary with respect to the normalisation of 
the zeroth order continuum strong residuals are present below 
10 keV. Keeping and Ai^ tied to each other, we found that 
a cross-normalization factor of ~ 3 is indeed required to 
reproduce the shape of the 4-10 keV continuum and the inten- 
sity of the Fe K emission lines, where the reflected component 
dominates. This forces the inclination angle between the axis of 
the torus and the observer's line of sight to a grazing value (~ 60°). 

Although we cannot rule out this scenario, we must allow 
for the possibility of a different geometry taking into account 
all the information that we obtained from the historical X-ray 
observations of NGC 4507 which suggested that: a) the column 
density of the line of sight (los) absorber varies; b) the flux of the 
Fe Kof remains rather stable, suggesting the presence of a constant 
and distant reprocessor and c) the primary continuum is also 
variable. Physically, the situation we want to model corresponds to 
a patchy reprocessor in which the reflected continuum is observed 
from reflection in matter on the far-side of the X-ray source, 
without intercepting any other "clouds," and the zeroth-order 
continuum corresponds to extinction by clouds in the line-of-sight. 
In practice this corresponds to allowing the column densities of 
the zeroth-order and reflected continua to be independent of each 
other; we thus followed the methodology discussed by lYagoobl 
I2OI2I applied to the modelling of the broadband X-ray emission 
of NGC 4945. By decoupling these two components, we can 
also allow the reflected and transmitted components to have a 
different temporal behaviour We can do this by decoupling the 
inclination angle parameters for the line-of-sight (zeroth-order) 
continuum passing through the reprocessor and for the reflected 
continuum from the reprocessor and allowing the column densities 
responsible for the reprocessing of the primary emission to be 
independent. The reflected continuum (and the fluorescent line 
emission, which is tied to it) is not extinguished by another column. 

The inclination angle of the zeroth-order component is now 
irrelevant (so it is fixed at 90 degrees), and the inclination angle 
for the reflected continuum is fixed at degrees because the effect 



of the inclination angle on the shape of the reflected continuum is 
not sufficiently large (in terms of spectral fitting) if the reflected 
continuum is observed in reflection only. Furthermore, since we are 
trying to model a patchy reprocessor as suggested from the column 
density variations, the inclination angle may not be meaningful. 
We note also that we cannot interpret the ratio between the normal- 
isations of the reflected and transmitted components simply as a 
covering factor of the reprocessors. Although the constant in front 
of the reflected continuum (j4r) does contain some information 
on the covering factor, that information cannot be decoupled from 
the effect of time delays between variability of the direct X-ray 
continuum and the reprocessed X-ray continuum. This is because 
the light-crossing time of the reprocessor is likely to be much 
longer than the direct X-ray continuum variability timescale, so 
the magnitude of the reflected continuum corresponds to the repro- 
cessed direct continuum that is averaged over a timescale that is 
longer than the reprocessor light-crossing time. This decoupling of 
the MYTORUS model is close to the standard procedure used while 
fitting with PEXRAV plus an absorbed power-law component. 
However there are several differences, in particular the column 
density of the reflector is also a free parameter and the Fe emission 
line intensities are calculated self consistently. 

The model then yields a x^/d.o.f. = 392.5/359 and a mean 
line-of-sight column density of iVn = 9.4to J x 10^^ cm"^, while 
the angle- averaged column density of the reflector (out of the 
line-of-sight) is Nh = 2£to%'X-lQ'^''^ cm"^, where also the Fe 
K emission lines are produced. The photon index is found to be 
r = I.68I0 Q3 normalisation of the primary continuum is 

l.gilo 25 X 10"^ ph cm"^ s"\ We also allowed the constant for 
the normalisation of the reflected continuum (Ar) to vary and we 
found that it is consistent withl(ylR = l.l±0.2). Finally, we note 
that now the measured velocity broadening of the Fe Ka emission 
line is a„ < 29 eV. 

4.1.2 Mytorus model for Suzaku, XMM-Newton and BeppoSAX 

We then applied the same model to the XMM-Newton observation. 
For simplicity, since there is no evidence of variability of the soft 
X-ray emission and taking into account that the Gaussian emission 
lines plus the thermal components are simple phenomenological 
models, we decided to keep fixed the main parameters of the latter 
to the Suzaku best-fit model. Furthermore, since we lack of simul- 
taneous observation above 10 keV we also fixed the photon index 
to the one measured with Suzaku and for simplicity at first we kept 
the constant of the relative emission line component (Ai^) fixed to 
1. We found that the out of los column density was comparable 
to the one measured during the Suzaku observation (TVh = 3.5 
_Q 4XlO cm ) while the los absorbing column density was 
A^H = 4.710.3 X 10^^ cm-2(xVdof=600.5/498). 

In contrast with the previous modelling with PEXRAV, we 
can now attempt to investigate also the relative intensity of the 
zeroth-order and reflected components. We found that the intensity 
of the zeroth order changed from (l-Olto^s) x 10"^ ph cm"^ s"^ 
to (1.56l|5:?2) X 10~^ ph cm"^ s"^ during the XMM-Newton and 
Suzaku pointing respectively. This suggests that no strong variation 
of the primary continuum is required to explain the observed 2-10 
keV spectral differences, but the main driver of the variations is 
the change in the column density of the line-of-sight absorber 

Finally, we applied the same model to the 3 BeppoSAX ob- 
servations, allowing also the photon index to vary, and we found 
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Table 3. Summary of the X-ray emission lines detected in tlie 6-8 keV energy range. Tlie energies of the lines are quoted in the rest frame. Fluxes and 
identifications are reported in column 2 and 3. The EW are reported in column 4 and they are calculated against the total observed continuum at their 
respective energies. In column 5 the improvement of fit is shown with respect to the continuum model, the value for the model with no lines is x^/d.o.f. 
=1918.9/375. 



Energy 
(keV) 
(1) 


Flux 

(10~''phcm-2 s-i) 
(2) 


ID 

(3) 


EW 
(eV) 
(4) 


Ax' 
(5) 


b.4U8_(j 


52.4«;I 


FeKa 


490l« 


1427.4 


'■'^'-0.02 




FeK^ 


8llf3 


82.0 


6 73 +°-''^ 
-0.06 


•'■^-1.6 


Fe XXV 


97 +15 
^' -13 


8.1 


'•-'"-0.08 


^■-'-1.4 


NiKa 




6.7 



Table 4. Comparison between the best fit values for the continuum and Fe Ka emission line for the Suzaku XMM-Newton and BeppoSAX observations. The 
fluxes are corrected only for Galactic absorption. 



Parameter 








Suzaku 


XMM-Newton 


SAXl 


SAX2 


SAX3 


DATE 








2007-12 


2001-01 


1997-07 


1998-07 


1999-01 


A^H (10^^ cm 










5 0+0'2 

^■'J-o.s 


7 0+0-4 
'-"-0.5 


6 2+0-7 


7 9+2.9 
'-^-1.6 


THard 








^ O-I+0.04 

l-»l-0.04 


1 76+0-05 
1^- 'D-0.04 


, 77+0.03 
l-"-0.05 


, 79+0.07 
'-'^-0.1 


1 6+°-i 


Normalisation 


(10 


^ ph cm ^ s- 


-1) 


64+° " 
'^■°*-o.ii 


-, r,c:+0.06 
^■^'-'-0.06 


9 ,9+0.22 
^-^^-0.33 


1 70+0-19 
^-'"-0.23 


62+0-48 

"■"^-0.23 


Tsoft 








o O+0.2 
^■*-0.2 


T 7+0.1 
^-'-0.1 








Normalisation 


(10 


""^ ph cm^'^ s" 


-1) 


1- 'U_o.09 


-, 07+0. 06 
-0.05 








^pcxrav (10 


2ph 


cm~^ s~^) 




1 1Q+0.09 
^■^^-0.09 




1 6+0-2 


1 9+0.3 
1-^-0.3 


1 0+0-5 


Fe Ka (keV) 








fi /ins+0 005 
b.40S_(, 


6.39lO-«^ 


6 39+0.10 
"--''^-0.09 


SR+0-11 


6 42+0-07 
°-^^-0.07 


/PE (10-5ph 


cm~ 


2 s-1) 




5 24+0-51 
^■^^-0.37 


4 36+0-43 


4 7+1.7 
^-'-1.8 


6 7+2-0 

°-'-2.2 


5 9+1-5 

■'-^-1.5 


EVKpE (eV) 








490+40 


190+20 
^^"-100 


140+80 

'^"-60 


99C+65 
^^J-90 


400+90 
^"'^-120 


F (0.5-2 kcV) 


(10 


— 1 ^ —9 

^ erg cm 




~ 4.8 


~ 4.9 








F(2-10 kcV) (10 


n —2 
erg cm s 




~ 0.6 


^ 1.2 


~ 1.8 


~ 1.6 


~ 0.8 


L(2-io kcV) (10'*^ erg s ^) 




~ 0.9 


~ 1.5 


~ 2.8 


^ 2.2 


~ 1.0 



that the column density of the out of the los absorber remained sta- 
ble and it was comparable (within the errors) to the one measured 
with Suzaku and XMM-Newton, (Nn = 2.6 +°-^xl0^3 cm"^ 
Nh = 2.4 15j-|xl0^3 cm-2 ^^^j j^^ ^ 2.4 tJ-^xlO^^ cm-2 foj. 
SAXl, SAX2 and SAX3 respectively). The column density of the 
los absorber varied with a similar trend as the one measured with 
the PEXRAV-based model (Nh = 6.3 x 10^^ cm"2^ ^ 55 



+0.5 , 



; +0.1 , 



_o.5 X 10^'' cm^^and Nh = 8.6 ; 
and SAX3 respectively), the intensity of the primary continuum 
components also varied. In particular, the intensity of the primary 
continuum was higher in the SAXl and SAX2 observation and in 
SAX3. However, due to the lower statistics of the BeppoSAX data 
these measurements have a large eiTor which prevent us from deriv- 
ing a clear picture. We note that there is no evidence for a variation 
of the photon index. 



5 DISCUSSION AND CONCLUSIONS 

Detailed X-ray spectral analysis of the Suzaku data confirms the 
complexity of the X-ray emission from NGC 4507. Thanks to 
the wide-band spectrum covering from 0.6 keV to 70 keV, we 
have now obtained the most reliable deconvolution of all the 
spectral components. The X-ray continuum is composed of three 
components: a heavily absorbed power-law component, a reflected 
component and a weak soft scattered component. We analysed the 
Suzaku and the historical observations of NGC 4507 adopting both 



a standard model as well as a new toroidal reprocessor model. With 
either of these two approaches we found that during the Suzaku 
observation the 2-10 keV emission of NGC 4507 was dominated 
by the reflected emission, while above 10 keV the spectrum was 
dominated by the highly absorbed transmitted component. 



The soft X-ray emission can be well described with a 
superposition of a power-law component, which is considered to 
be the scattered light from ionised, optically-thin gas, and several 
emission lines. These e mission lines were al ready detected with 
the ASCA observation ( IComastrietal.lll998l). As alre ady shown 
with the XMM-Newton observation jMatt et al.l l2004h . the wide 
range of ionisation implied by these lines is indicative of the 
presence of at least two photo- or collisionally- ionised emitters. 
This is confirmed by a recent deep XMM-Newton observation 
obtained by our group within a monitoring program of NGC 4507. 
The analysis of the single observations showed also that as seen 
in other Seyfert 2s the soft X-ray emission did not vary, implying 
that the emitters responsible for this emission are loc ated outside 
the variable X-ray absorber jMarinucci et al. |[2012bl . Wang et al 
private communication). 
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Table 5. Comparison between the best fit values for the continuum and A^h for the Suzaku, XMM-Newton and BeppoSAX observations when fitted with the 
MYTORUS model. The fluxes are corrected only for Galactic absorption, while the luminosities are corrected also for the intrinsic absorption, which includes 
the effect of the Compton-down scattering. 



Parameter 


Suzaku 


XMM-Newton 


SAXl 


SAX2 


SAX3 




2007-12 


2001-01 


1997-07 


1998-07 


1999-01 


r 

Normahsation (10~^ ph cm~^ s— 1) 

N-a transmitted (10^^ Cm-^) 
A^H reflected (10^^ Cm'^) 


1 Q1+0.16 

q OO+0.25 
».oo_Q 24 

Z.OO_Q 47 


1.68^' 
1 56+0-°5 

^.OD_g 27 

o ,54+1.26 


1 62+''°^ 

^■''^-0.35 
.33+0.29 

2 64+0.88 


1 ™_0.05 
9 qo+0.49 

"■™_0.46 

2 4Q+0.85 


^■"-'-0.07 

^■^^-0.61 
,0+0.10 

+ 1.29 
^■•"-1.05 


F(2-10 keV) (10"" ergcm-2 s"!) 


~ 0.6 


~ 1.2 


~ 1.7 


~ 1.6 


- 0.9 


L(2-iokeV) (10«ergs-l) 


~ 2.5 


~ 2.1 


~ 4.0 


~ 4.1 


~ 3.0 



5.1 The X-ray absorber 

In the last two decades NGC 4507, which is one of the X- 
ray brightest and nearby Seyfert 2 galaxies, has been ob- 
sei^ed several times with all the different X-ray observatories; 
NGC 4507 displayed an observed 2-10 keV flux ranging from 
0.6 — 1.3 X lO^i^erg cm^^ s^^; furthermore these X-ray ob- 
servations showed long-term A^h variability, which changes by a 
facto r of 2, and also possible variability of the intrinsic continuum 
jRisa liti 2002). The Suzaku observation caught the source with 
a low observed 2-10 keV flux, similar to the last BeppoSAX 
observation (SAX3), which was also characterised by the highest 
measured column density of the X-ray absorber (8 — 9 x 10^'' 
cm~^). Our analysis of the XMM-Newton BeppoSAX, and Suzaku 
observations shows that, independently from the assumed model 
(i.e. MYTORUS or the standard PEXRAV one) the X-ray absorber in 
NGC 4507 varies from ~ 5 x 10^^ cm~^ to ~ 9 x 10^^ cm"^ on 
~ 6 months time scale between the observations. 

Assuming a spherical geometry for the obscuring clouds, and 
that they ar e moving with Kep lerian velocities (as in the case of 
NGC 1365: iRisaliti et alj|2007h , then there is a first order simple 
relation between the crossing time of such obscuring cloud, the 
linear dimension and di stance of the obscuring c loud and the size 
of the X-ray source (see lMarinucci et ani2012bh . From this work, 
since no variability is found on a short time scales, during the long 
Suzaku observation, we can infer that the variability occurs on a 
time scale between 2 days (as the elapsed time of the Suzaku obser- 
vation was ~ 180 ks) and six months (as the elapsed time between 
the second and the third BeppoSAX observa tions). Following th e 
same argument proposed for NGC 1365 jRisaliti et all l2007h . 
where similar assumptions are used, then a possible scenario could 
be that the variable absorber is located a distance greater than 
0.01 pc from the X-ray source (i.e. not closer than the Broad Line 
Region). Thus a possible scenario, where the classical uniform 
absorber still exists is that there are multiple absorbers. One is the 
classical and uniform pc-scale absorber (i.e. the torus), which is 
responsible for the Fe K emission line and the constant reflected 
component; while the A^h variability requires the presence of a 
second and clumpy absorber that could be coincident with the 
outer BLRs. We note that the observations presented here do 
not span the possible time-scale expected for a variable absorber 
located between the BLR and the pc-scale absorber. 

However a simpler scenario could be that the variability 
is due to a certain degree of dump iness of a single absorber 
itself, as for the clumpy t orus model jNenkova et al.|[20o3 . |2008| ; 
lElitzur & Shlosmanlliooeh . In this scenario, if the distribution of 



clouds remains rather constant, we can have a constant reflected 
component while at the same time changes of the line-of-sight 
TVh. Recently our group obtained an XMM-Newton monitoring 
campaign of NGC 4507 consisting of six observations performed 
every 10 days, this monitoring confirmed that the A'h variabil- 
ity occurs on relatively long time scale (between 1.5 and 4 
months), supporting the scenario of a clumpy pc-scale absorber 
jMarinucci et al. ir2012bl) . A long term variability or lack of certain 
components cannot be simply used to infer the location of the 
reprocessor itself in terms of time delays between the reflected 
component and the primary continuum. As discussed in several 
works on the s pectral variability of well m onitored Seyfert Is (e.g 
MCG-6-30-15 iMiller et alj|2008l ; Mrk 766 lMiller etallbOOTl) . the 
apparent variability of the observed continuum level could be due 
to changes of the los covering fraction and not of the intrinsic 
continuum level. Thus the reflected emission would remain rather 
constant even for closer in and variable absorbers. 

Stronger evidence for the presence of a pc-scale reflector can 
be derived from the analysis of the Fe Ko? emission line profile. 
We note that Fe Ka emission line is rather constant and narrow 
with no evidence of an additional strong broad component; its 
width, measured with Suzaku, is cr < 30 eV (or av < 1400 km 
s^i) and consistent with the Chandra upper limits, (t„ < 1200 km 
s^i {FWHM < 30 00 km s~i). The m easured FWHM of the 
Ha (~ 5000 km s~i; lMoran etai]|2000l) is marginally higher than 
the Chandra upper limit on the Fe Ka FWHM. This suggests 
that the Fe Ka is produced either in the outer part of the BLRs 
or in a pc-scale absorber; in agreement with a scenario where 
there is a distant and stable reprocessor, which could be identified 
with the c lassic al torus. A similar result has been presented by 
IShu et alj lioiH) , where the authors compared the FWHM of the 
Fe Ka emission line (measured with the Chandra HETG) of a 
sample of Seyfert 2s (including NGC 4507) with the FWHM 
of the optical lines. They suggested that the Fe Ka emitter is a 
factor 0.7-1 1 times larger than the optical line-emitting region and 
located at a distance of about r ~ 3 x 10* rg (where rg is the 
gravitational radius defined as rg = GM/c?). In particular from 
the e stimated black ho le mass for NGC 4507 of Mbh ~ 4 x lO'^ 
M0 jSian & Gill2007l) and assuming Keplerian motion, the upper 
limit on its FWHM implies that the Fe Ka emission line is 
produced at a distance r > 0.02 pc from the central BH. We 
note that assumi ng the large r BH mass (Mbh ~ 2.5 x 10^ 
Mq) reported by I Winter et "all l l2009l) would place the absorber at 
a distance R > 0.06 pc and the location of the Fe Ka at r > 0. 1 pc. 

In terms of the global picture for the location and structure of 
the X-ray absorbers, we have now several examples of obscured 
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AGNs with short-term variation unveiling that a significant 
fraction of the absorbing clouds are located within the BLR. 
However there is also evidence for the presence of a pc-scale 
absorber as predicted in the Unified Model of AGNs. This absorber 
is confirmed by the ubiquitous presence of the narrow Fe Ka 
emission line and the Compton reflection component, which 
do not show strong variability between observations even in 
case of a variable intrinsi c continuum and/or v a riable neutral 
absor ber (e.g. N GC 7582 IPiconcelli et aP l2007l: I Bianchi et al 



2009L NGC 4945: lMarinucci et al.ll2012al : lYaaoobll2012l : lltoh et al 
20081) . Another piece of evidence for the presence of a distant 



reprocessor comes from the comparison between the Chandra, 
XMM-Newton and Suzaku observations of the brig ht Seyfert 
2 NGC 4945 l lMarinucci et"ai] 1201 2al : lYaooobI 12013) . where a 
detailed spectral, variability and imaging analysis unveiled that the 
emitting region responsible for the Fe Ka line and the Compton- 
scattered continuum has a low covering factor and it is most likely 
located at a distance > 30 — 50 pc. In this framework the relatively 
long-term absorption variability shown by NGC 4507 confirms 
that the location and structure of the X-ray absorber is complex 
and that absorption in type 2 AGNs could occur on different scales 
and that there may not be a universal single and uniform absorber. 

Finally, by adopting the standard PEXRAV model, even with 
the broad band X-ray observations available for NGC 4507, we 
cannot assess the role of the possible variability of the primary 
continuum or estimate the column density of the reprocessor 
responsible for the Fe Ka emission line. Interestingly, the fit with 
the decoupled MYTORUS model, which can mimic either a patchy 
toroidal reprocessor or a situation where there are two reprocessors 
(one seen in "transmission", dominating the high-energy spectrum, 
and one seen in reflection), allows us to measure these column 
densities and the possible variability of the primary continuum. 
Although the decoupled MYTORUS model closely resembles the 
classical combination of PEXRAV (slab reflection component) 
and an absorbed power-law component, the column densities 
of both the reprocessors ("reflector" and "absorber") are treated 
independently and self-consistently with the emission of the Fe 
K line. Table |5] shows that the line-of-sight obscuration of the 
reprocessor seen in transmission varies by AA^h ~ 5 x lO'^^ 
cm^'^, while there is a constant reprocessor with a column density 
of ~ 2 — 3 X lO^'^cm"^, which is the one responsible for the 
Fe Ka emission line. The intrinsic X-ray luminosity ranges from 
i{2-iokcV) ~ 2.1 X 10*=* erg s'^o L(2_io kcV) ~ 4.1 x 10*^ 
erg s^ ^ . While we observe intrinsic variation of primary power-law 
intensity, the A'h variablity drives the spectral changes between 
2-10 keV. Conversely, the reflection and emission line components 
are not observed to vary. 

We note that the behaviour of the line-of-sight column and 
reflection fraction with respect to the intrinsic continuum going 
from the BeppoSAX to the Suzaku observation strongly depends 
on the adopted model for the reprocessor. Indeed as can be seen 
comparing Table|4]and [5] by adopting the combination of PEXRAV 
and an absorbed power-law component, we would infer that 
NGC 4507 was intrinsically brighter and less obscured during the 
XMM-Newton observation. However, no such trend is inferred 
with the MYTORUS model, where the opposite is the case (i.e. the 
more absorbed Suzaku observation has the higher primary power- 
law normalisation). Only future monitoring campaigns with broad 
band observatories such as ASTRO-H (i.e. with instrument with an 
high effective area above 10 keV as well as high spectral resolution 



at the Fe Ka line) will allow monitoring of sources like NGC 4507. 
These observations will allow to investigate the variability of the 
harder continuum simultaneously providing a detailed investiga- 
tion of the profile of the Fe Ka emission line, thus establishing 
the geometry and location of the "stable" and variable reprocessors. 
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